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基于 投影 法 求解 不 可 压缩 流 的 高 精度 紧 致 格式 * 


T 坪 ， 田 振 大 t 
(复旦 大 学 力学 与 工程 科学 系 ，. 上 海 200433) 
摘 E: 本 文 建立 了 一 种 基于 投影 法 的 求解 不 可 床 缩 Navier-Stokes (N-S) 方程 的 高 精度 紧 致 差分 格式 。 
该 方法 时 间 上 上 采用 Kim 和 Moin 二 阶 投影 法 离散 ， 空 间 上 采用 高 精度 紧 致 格式 离散 ， 并 提出 了 一 
种 新 的 离散 压力 边界 的 紧 臻 格式， 同时 对 计算 结果 进行 分 析 以 验 让 该 投影 法 的 精度 和 格式 稳定 
性 。 文 中 Taylor 涡 列 数值 计算 结果 表明 ，Kim Moin 投影 法 能 使 得 讨 力 场 和 速度 场 均 达到 时 间 
一 阶 精度 ， 且 高 精度 紧 致 格式 投影 法 也 具有 空间 高 阶 精 度 。 了 驱动 方 腔 数值 模拟 结果 显示 ， 本 文 
对 N-S 方 程 的 离散 格式 具有 很 好 的 可 靠 性 ， 适 用 于 对 复杂 流体 流动 的 小 尺度 问题 的 数值 模拟 和 研 


究 。 
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1 引言 
描述 不 可 压缩 非 定常 流体 的 控制 方程 为 Navier-Stokes (N-S) 方程 组 ， 其 无 量 纲 形式 可 写 为 
O07 vv Eee 
a, tu Vi= VP+ ReY a, (1) 


其 中 如 为 速度 失 量 ，p 为 压力 ，Re 为 雷诺 数 。 

N-S 方 程 组 不 是 完全 时 间 发 展 型 方程 ， 在 每 个 时 间 步 上 速度 必须 满足 散 度 为 零 的 约束 条 
件 (2)， 并 出 现 压力 -速度 看 合 现 象 。N-S 方 程 组 的 求解 难点 在 于 如 何 将 速度 和 压力 解 艳 求解 。 
投影 法 的 基本 思想 是 通过 引入 一 个 无 需 满足 散 度 为 零 约束 条 件 的 中 间 变 量 对 压力 速度 解 看 
并 求解 ， 通 过 这 种 算 子 分 裂 后 分 步 求 解 的 方式 能 够 高 效率 地 计算 非 定 常 问题 ， 因 而 得 到 了 
非常 广泛 的 应 用 。 域 初 由 Chorin! 提出 的 投影 法 只 具有 时 间 一 阶 精度 ， 之 后 很 多 学 者 发 展 了 
提高 时 间 精 度 的 投影 法 。Kim 和 Moin 提 出 了 一 种 具有 了 时间 二 阶 精度 的 无 压力 投影 法 外 ，Van 
kan 和 Bel 等 人 提出 了 一 种 简化 中 间 变 量 边 办 条件 的 压力 投影 法 84， 刘 森 儿 ， 任 下 新 等 人 对 
投影 方法 研究 进展 进行 了 绿 述 ， 并 发 展 了 二 阶 迭 代 投 影 法 以 及 时 间 上 有 具 有 二 阶 精度 的 投影 
法 下 ?1，Brown，Cortez 等 人 对 不 同 的 投影 法 作 了 有 具体 的 精度 分 析 并 提出 了 改进 的 方法 加 。 

在 流体 模拟 中 ， 很 多 问题 都 是 小 尺度 流动 。 相 对 低 精 度 格 式 ， 高 阶 紧 致 格式 由 于 具有 精度 
高 的 优点 而 能 吏 好 地 模拟 小 尺度 问题 名。 对 于 光滑 问题 ， 紧 致 格式 还 能 降低 对 网 格 基 架 点 的 要 
求 ， 从 而 提高 计算 效率 。 相 对 于 对 称 格式 ， 迎 风格 式 具 有 的 耗 散 性 可 以 压制 数值 解 中 非 物 理 
高 频 振荡 ， 因 此 迎风 格式 得 到 很 大 的 发 展 ， 并 出 现 了 各 种 高 精度 迎风 格式 10 坪 。 为 了 提高 中 
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近 解 率 ， 我 们 将 五 阶 紧 致 迎风 格式 和 四 阶 紧 致 迎风 格式 进行 加 权 得 到 了 一 种 组 合 迎 风 紧 致 格 
AN2, Abide 和 Viazzo 以 Kim 和 Moin 无 压力 投影 法 为 基础 发 展 了 一 种 高 精度 对 称 紧 致 格式 的 
投影 法 HB， 计算 结果 显示 压力 场 和 速度 场 在 空间 上 能 达到 高 阶 精度 ， 但 在 时 间 层 上 压力 场 却 
只 有 一 阶 精度 ， 但 Brown!) 通过 正则 模 态 分 析 认为 该 投影 法 的 压力 场 具 有 时 间 二 阶 精 度 。 本 文 
基于 无 压力 投影 法 四， 发 展 了 一 种 不 可 压缩 流动 N-S 方程 组 的 高 精度 紧 致 差分 格式 。 该 方法 在 
时 间 上 采用 Kim 和 Moin 投影 法 离散 ， 空 间 上 采用 紧 致 格式 离散 ， 并 提出 了 一 种 新 的 压力 边界 
离散 格式 。 本 文 所 采用 的 数值 算 例 中 ，Tayior 涡 列 数值 计算 结果 很 好 地 验证 了 Kim 和 Moin 投 
影 法 能 使 得 压力 和 速度 在 时 间 上 均 达 到 一 阶 精 度 ， 且 在 空间 上 达到 了 四 阶 精度 ， 驱 动 方 腔 算 例 
的 计算 结果 显示 本 文 格式 具有 很 好 的 计算 稳定 性 和 可 靠 性 ， 同 时 对 于 复杂 流体 流动 的 小 尺度 问 
题 也 有 很 好 的 模拟 结果 


2 ”控制 方程 的 离散 和 数值 算法 


2.1 控制 方程 的 时 间 离 散 
不 可 压 N-S 方 程 在 时 间 上 的 二 阶 精度 离散 格式 为 


antl _ gn + (@- Va)rt? = -Vpt + meV + a), 
At (3) 
V-t =0, 
其 中 人 t 为 时 间 步 长 ， 对 流 项 (了. Vat: 可 采用 Adams-Bashforth 的 一 阶 格式 处 理 如 下 : 
(a-vat? = (a va- 3(U Vay, (4) 
速度 相应 的 边界 条 件 为 Dirichlet 边界 条 件 
TH Fes = apt, (5) 


投影 法 是 根据 Helmholtz-Hodge 分 解 定 理 ， 引 入 一 个 中 间 速 度 场 和 压力 场 ， 对 速度 和 
压力 解 耦 后 分 步 计算 ， 从 而 把 (3) 式 的 求解 分 为 预算 步 ， 修 正 步 以 及 最 终 步 。 本 文采 用 的 
是 Kim 和 Moin 提出 的 二 阶 投影 法 ， 具 体 步骤 如 下 。 

步骤 1 预算 步 。 引入 中 间 速 度 场 并 求解 

这 一 人 
At 


1 


21% =n 
= zV (T +T"), (6) 


十 二 (a. va” — 5 (i vi! = 


nv). =n grt 
lan len? (7) 


tee T Jan i Ui + At: Ve" log 


FE az, A Su yf ARREA AU) FR Jy Hp Bt HERED BE BUE Dy EY A EA 
R (7) Æ Kim Al Moin 提出 的 中 间 速 度 边界 条 件 。 
步 又 2 ”修正 步 。 通 过 下 式 对 预算 步 得 到 碎 的 进行 修正 


Cat AL vet, (8) 


V-at! =0, (9) 
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其 中 Gm+1 ERA prt 的 某 种 近似 ， 我 们 称 为 中 间 压 力 场 ， 应 满足 


At 


nn 十 村 = nt+1_ 7 
K ? Re 


wager. (10) 


修正 后 速度 we 的 散 度 为 零 ， 既 满足 连续 性 方程 ， 且 满 足 动量 方程 。 
对 式 (8) 两 端 取 散 度 ， 并 结合 (9) 得 到 压力 Poisson 方程 ， 并 添加 边界 条 件 ， 得 到 


ygn = E E, 
agt! (11) 
ðn laa r 
步骤 3 Be RUTAN A, AAP t WR peta 
atl = g At VE, 
12) 
n+} 一 o” ,_ Ot 2 pm 十 1 ( 
prt? = gn mv gntt, 
2.2 ”控制 方程 的 空间 离散 
控制 方程 中 的 对 流 项 我 们 可 采用 迎风 格式 来 壳 近 ， 对 流 项 可 用 如 下 的 分 裂 形式 表 不 
Ou edt , Ou 
‘ae ae On’ (13) 
a lul juj 
"WE u+ lu —_ u-u 
u = 2 u = 2 ’ 
iB ut 2, use 4 Rut Ft, u FO ima, mh 表 未 一 阶 偏 导 数 Se 在 点 ;处 的 通 近 值 ， 
我 们 采用 如 下 的 四 阶 迎 风格 式 通 近 
1 7 1 
Ft, +2FF = 7(- guint 二 i 一 suit ) + hsi, 
(14) 


= = 171 7 
2F; + Fip = 5 (gu — 4ui + suit ) — hSi, 


其 中 8 KROMA S SS CEG ANAL, RAN Pade fk SLI» 74 在 计算 中 
涉及 到 的 不力 梯度 项 的 求解 ， 我 们 采用 一 阶 导数 四 阶 对 称 Padé FRECHE SCRE 
3 ” ”Poisson 方程 的 数值 算法 


在 投影 法 中 ， 最 关键 的 预算 步 和 修正 步 都 可 以 归纳 于 求解 Poisson 方程 。 动 量 方程 (6) 和 压 
力 Poisson 方 程 (11) 可 以 统一 写成 如 下 Poisson 方程 形式 


ve =f, (15) 
barih ROS RRHRWE ©), T 的 边界 条 件 为 Dirichlet WHA. AINE 


© x = 3) An 1. —\n—1 2RE mn 2m 
Pi S = 2Re [5@-va g(t We) | age Vu", 
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二 $ 时 ， 式 (15) 表 示 动 量 方 程 (11)，9 的 边界 条 件 为 Neumann 边界 条 件 。 相 应 地 
V g* 
人 


3.1 Poisson 方程 的 高 精度 紧 致 差分 格式 
为 了 提高 计算 效率 以 及 保证 压力 和 速度 的 精度 ， 我 们 采用 四 阶 九 点 紧 致 格式 H4 对 方 
程 (15) 进行 离散 ， 其 四 阶 九 点 紧 致 离散 格式 如 下 


4 8 4 
4924.4 308, — 206 TEA ODEL] (16) 
t=1 i=5 i=1 


Hosa, f= 5+ TETRA (OWE 


AEN 8 Re 
hn) + 7 — (20+ 8h? Jag = E (XS +850), (17) 
#4 & Dirichlet 边界 条 件 ， 采 用 ADI 2S ( 松 驰 ) 迭代 法 对 式 (17) 进行 求解 。 
3.2 TEJ Poisson 方程 的 边界 处 理 
由 于 压力 Poisson 方 程 (11) 具有 Neumann 边 界 条 件 ， 为 了 保证 格式 的 整体 精度 以 及 稳定 
性 ， 我 们 对 Neumann 边界 条 件 OS on = 0 构造 了 一 种 新 的 紧 致 格式 。 例 如 ， 在 右边 界 及 其 
右上 角 点 的 离散 格式 为 
1 h? 3 
bi41,5 + 5 (Gi,541 + ij-1) — 261,53 = > f +O?) 
2 a 2 (18) 
Pi+1,j + Pig—-1 — 261,75 = z’ + O(h?), 


A (18) 是 一 种 具有 三 阶 精度 的 紧 致 离散 格式 ， 在 后 文 的 数值 计算 中 也 显示 ， 该 离散 格式 在 边界 
点 并 无 损失 精度 从 而 保证 了 算法 的 整体 精度 。 结 合 内 点 的 离散 方程 (16) 即 得 到 一 封闭 的 方程 
4, RH ADI 2% ( 松 驰 ) 迭代 法 进行 求解 。 

3.3 压力 Poisson 方 程 解 的 唯一 性 

对 于 具有 Neumann 边界 条 件 的 压力 Poisson 方程 


V6 =f, 


3$ 
On lan 


其 解 并 不 唯一 。 采 用 文献 [15] 的 策略 ， 我 们 可 以 添加 一 个 附加 条 件 ， 使 方程 (19) 变 为 


(19) 


一 b 


vV = 也 


Ob) 
én an — 0, (20) 


2 9i,; = 0, 
增加 附加 条 件 $i,; = 0 后 的 方程 (20) RAE fk, UAT EE ae E Je A BR (1905 
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4 RAAH 


4.1 ”Taylor 涡 列 问 题 
Taylor 涡 列 问题 有 如 下 精确 解 


u(z,y,t) = — exp(—2tr? /Re) cos(nx) sin(ry), 
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v(x, y, t) = exp(—2tn?/Re) sin(rz) cos(zy), (21) 


plz, y,t) = -5 exp(—4tr?/Re) [ cos? (rz) + cos? (ny)], 


该 问题 经 常 被 用 来 判断 数值 方法 精度 和 算法 的 可 靠 性 。 选 取 空 间 区 域 为 0 < zx, y <1, Re = 
10。 为 了 判断 本 文 格式 的 时 间 精 度 ， 我 们 选择 不 同 的 时 间 步 长 来 计算 在 空间 步 长 为 h = 
1/58, t = 0.3 时 刻下 的 最 大 绝对 误差 。 表 1 结果 显示 ， 本 文 的 投影 法 格式 使 速度 场 和 压力 场 均 
达到 了 二 阶 精度 ， 压 力 场 的 收敛 精度 略 小 于 速度 场 的 精度 。 同 时 ， 为 了 判断 本 文 格式 的 空间 精 
度 ， 我 们 选取 时 间 步 长 为 + = 1/2000， 在 t= 0.3 时 刻 不 同 空间 步 长 下 的 最 大 绝对 误差 。 表 2 绪 


果 显 示 ， 本 文 格式 能 够 保证 压力 场 和 速度 场 均 能 达到 空间 四 阶 精度 。 


Rl: ”空间 步 长 hh 二 1/58, t= 0.3 时 刻 的 最 大 绝对 误差 


v rate 
|2150 | 了 


了 rate 
2.1210E-4 
1.5159E-4 1.842 
8.4082E-5 | 2.049 
5.3771E-5 2.003 
3.7779E-5 1.936 


表 2: 时 间 步 长 r = 1/2000, t= 0.3 时 刻 的 最 大 绝对 误差 


一 ”| 
Ce | sema semma | | | 
Co | mume | as% | 275006 | same | L005.4 | 
Co | eoseae-7 | 3770 | 600087 | 3773 | 210565 | 


4.2 驱动 方 腔 问 题 
了 驳 动 方 腑 流 是 检验 数值 方法 稳定 性 的 常用 算 例 。 本 文采 用 文献 16] 中 的 修正 驱动 方 腔 ， 


1/58 


1/2000 


其 


ERRUR REKU = 16r?(1-— r), v= 0， 其 余 壁 处 的 水 平 速度 和 垂直 速度 均 为 0。 本 文 
取 时 间 步 长 为 + = 1/500， 空 间 网 格 为 128 x 128 进 行 计算 。 计 算 结 果 与 Shenl16] 的 结果 非常 相 
近 ， 对 主 涡 和 小 涡 的 大 小 及 位 置 都 有 非常 好 的 捕 提 。 表 3 为 本 文 格式 所 得 到 的 涡 心 大 小 及 位 


置 与 文献 上 16] 的 对 比 。 结 果 显 示 ， 本 文 对 主 涡 和 一 次 涡 的 物理 特征 有 很 好 的 模拟 结果 。 图 


1 Ail 


图 2 分 别 为 Re = 400, 1000, 2000 和 Re = 5000, 7500, 10000 时 方 腔 流动 达到 稳定 后 的 流 线 图 。 
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表 3: 主 涡 及 底部 二 次 涡 位 置 及 流 吕 数值 与 文献 的 比较 


原始 涡 心 处 (坐标 ) | 左下 角 涡 心 处 (坐标 ) 


AF A skh y (BFR) 
1.98E-4(0.922,0.094) 

2.18E-4(0.914,0.094) 
5.68E-4(0.922,0.094) 
6.41E-4(0.914,0.094) 
8.08E-4(0.922,0.094) 
9.26E-4(0.914,0.094) 
7.75E-4(0.922,0.094) 
9.09E-4(0.914,0.094) 
8.05E-4(0.914,0.094) 
3.25E-4(0.922,0.094) 
4.17E-4(0.914,0.094) 


400 | -0.0858(0.578,0.625) | 6.31E-6(0.031,0.047) 
1000 | -0.0871(0.547,0.578) | 8.28E-4(0.078,0.063) 
e -0.0873(0.539,0.570) 
2000 | -0.0877(0.531,0.547) | 3.54E-4(0.094,0.094) 
ea ee eee 
5000 | -0.0880(0.516,0.531) | 7.53E-4(0.094,0.094) 
B ren 


10000 | -0.0883(0.536,0.531) | 7.92E-4(0.094,0.094) 
-0.0893(0.516,0.531) | 7.94E-4(0.094,0.094) 


当前 


当前 
文献 [16] 


> 2 


图 2: Re=5000, 7500, 10000 方 腔 内 部 流 线 图 
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Æ XX} Navier-Stokes 方 程 在 时 间 上 采用 了 Kim 和 Moin 的 无 压力 投影 法 离散 ， 空 间 上 分 别 

采用 了 四 阶 迎 风 紧 致 格式 离散 。 在 Poisson 方 程 求解 离散 过 程 中 采用 了 九 点 四 阶 紧 致 格式 ， 并 
对 Neumann 边界 条 件 提出 了 具有 三 阶 精度 的 紧 致 差分 格式 。Taylor 涡 列 数值 结果 表明 ， 本 文 
所 提出 的 高 精度 紧 致 莽 分 投影 法 在 空间 上 ， 速 度 场 和 压力 场 均 达 到 了 四 阶 精度 。 在 时 间 层 上 ， 
压力 和 速度 在 整 场 中 也 都 达到 了 一 阶 精度 。 同 时 ， 驱 动 方 腔 的 计算 结果 也 显示 本 文 提出 的 格式 
具有 很 好 的 稳定 性 ， 对 复杂 流体 流动 有 着 很 好 的 数值 模拟 。 
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A Projection-based High-order Compact Scheme for Solving 
Incompressible Flows 


WANG Ping, TIAN Zhen-fut 
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Abstract: In this paper, based on the Kim and Moin projection technique, we develop a high-order 
compact difference method with spatial order 4 and temporal order 2 for the solution of the incom- 
pressible Navier-Stokes equations. In particular, we propose a new compact discretization formula for 
the Neumann boundary condition of pressure. Numerical experiments are performed to validate the 
accuracy and robustness of the presented scheme. The simulation results of the regularized driven 
cavity are in good agreement with the existing solutions in the literature. 
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